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ABSTRACT 

We present J' and K' imaging linear polarimetric adaptive optics observations of NGC 
1068 using MMT-Pol on the 6.5-m MMT. These observations allow us to study the 
torus from a magnetohydrodynamical (MHD) framework. In a Oi'S (30 pc) aperture 
at K', we find that polarisation arising from the passage of radiation from the in¬ 
ner edge of the torus through magnetically aligned dust grains in the clumps is the 
dominant polarisation mechanism, with an intrinsic polarisation of 7.0%±2.2%. This 
result yields a torus magnetic field strength in the range of 4—82 mC through para¬ 
magnetic alignment, and through the Chandrasekhar-Fermi method. The 

measured position angle (P.A.) of polarisation at K' is found to be similar to the P.A. 
of the obscuring dusty component at few parsec scales using infrared interferometric 
techniques. We show that the constant component of the magnetic field is responsible 
for the alignment of the dust grains, and aligned with the torus axis onto the plane 
of the sky. Adopting this magnetic field configuration and the physical conditions of 
the clumps in the MHD outflow wind model, we estimate a mass outflow rate ^0.17 
Mq yr“^ at 0.4 pc from the central engine for those clumps showing near-infrared 
dichroism. The models used were able to create the torus in a timescale of ^10^ yr 
with a rotational velocity of ^1228 km s“^ at 0.4 pc. We conclude that the evolution, 
morphology and kinematics of the torus in NGC 1068 can be explained within a MHD 
framework. 
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1 INTRODUCTION 


The torus is the cornerstone of the unified model | |Lawren^ 
1991 Antonuc'cl] 1993[ |Urry fc Padovani 19951 of Active 
Galactic Nuclei (AGNs). Specifically, this model posits 
that the observational differences between radio-quiet AGNs 
arise from an orientation effect. In this scheme, the AGN 
classihcation solely depends on the anisotropic obscuration 
of the central engine (black hole and accretion disc) by an 
optically and geometrically thick, dusty torus. Although the 
torus was postulated to be a “donut-shape” morphology, in 
this paper the term “torus” is used to denote the region sur¬ 
rounding the central engine where the obscuring material is 
located, with the precise morphology of that region still to 
be determined. 

Several ideas to explain the existence of the torus have 
been proposed. Some models explained the torus as an in¬ 
flow of gas from few tens of kpc. |Wada, Papadopoulos fc] 
Spaans (20091 presented numerical simulations of the inter¬ 


stellar medium to track the formation of molecular hydrogen 
forming an inhomogeneous thick disc around the central en¬ 
gine, identified as the torus. In a further 3D hydrodynamic 
simulation model, |Wada| ( |2012[ ) found that a vertical circu¬ 
lation of gas in the central parsecs drives the formation of 
an obscuring region along the equatorial plane. [Schartmann,| 
Krause & Burkert (20111 suggested the origin of the torus 


as dusty clouds falling from the host galaxy to the central 


engine. Gonversely, some models (jBlandford fc Payne |1982 

Krolik fc Begelman||19881 [Emmering, Blandford fc Shlos- 

man 

1992| Konigl fc Kartje||1994| [Kartje, Konigl fc Elitzur 

1999| 

Elitzur fc Shlosman|2006 

1 assume an outflowing wind 


confined and accelerated by the magnetic field generated in 
the accretion disc. In this scheme, the hydromagnetic wind 
can lift the plasma from the midplane of the accretion disc 
to form a geometrically thick distribution of dusty clouds 
surrounding the central engine. High angular resolution ob¬ 
servations suggest the infrared (IR) emitting regions of the 


torus to no more that a few parsecs in radius (i.e., Radom- 

ski et al. |j2003||Jaffe et al.||2004| 

Packham et al.||2005 Ma- 

son et al.j 20061 Radomski et al. 

20081 Ramos Almeida et 

al. 2009| [Raban et al. 20091 Ramos Almeida et al. 20111 

Alonso-Herrero et al.||20111 Burtscher et al. 

20131 Ichikawa 

et al. 2014 Markowitz, Krumpe fc Nikutta 

20141. Within 


these scales, the torus is in the region where the supermas- 
sive black hole (SMBH) and accretion disc activity dominate 
over the galaxy components. Thus, the torus origin is likely 
influenced by the central engine, which gives strong impetus 
to the outflowing wind model. 

Although large efforts have been made in the develop¬ 
ment of (magneto-)hydrodynamical (MHD) simulations of 
the outflow wind of AGNs, the magnetic field strengths in 
the torus are poorly constrained. The magnetic field can in¬ 
duce a preferential orientation of dust grains within the torus 
that can give rise to a measurable degree of polarisation (e.g. 
Davis fc Greenstein|195l]|Jones &: Spitzer|1967|[Jones, Klebe| 


fc Dickey|1992 1, with an expected linear degree of polarisa¬ 

tion <15% at near-IR (NIR: 1—5 fim) assuming standard 
disc models ( Efstathiou, McCall fc Hough|19^ l. Hence, IR 
polarimetric techniques give us a powerful tool to enhance 
the contrast of the polarised structures arising from aligned 
dust grains, helping us to constrain the magnetic fields, from 
those unpolarised components within the core of AGNs. Us¬ 


ing J, H and K„ imaging polarimetric observations on the 
3.9-m Anglo-Australian Telescope, |Lopez-Rodriguez et al.| 
( |2013[ ) estimated a magnetic field strength in the range of 
12—128 mG in the NIR emitting regions in the torus of IC 
5063. Through comparison with the magnetic field strength 
of ~2 mG estimated by the ratio of thermal and magnetic 
field pressure, they suggested that a magnetically dominant 
region is surrounding the central engine of IC 5063. They 
concluded that the magnetic field plays an important role 
in the torus evolution of IC 5063, favoring the MHD outflow 
wind model. 

NGC 1068 is an archetypal type 2 AGN, whose prox¬ 
imity (we adopt Ho = 73 km s“^ kpc“^, thus 1” = 60 pc) 
and high brightness make it an ideal target for polarimetry. 
The detection of polarised broad emission lines in the optical 


(0.35—0.70 fim) wavelengths of NGC 1068 by Antonucci & 


Miller (19851 was crucial to understand the AGN structure. 


Their detection is most readily interpreted through scatter¬ 
ing of the radiation from the central engine into our line 
of sight (LOS) by the ionisation cones. After correction for 
starlight dilution, the continuum polarisation was shown to 
be constant, ~16%, over the ultraviolet (Antonucci, Hurt & 


Miller 19941 and optical wavelengths (Miller fc Antonucci 


19931, with electron scattering being the dominant polarisa¬ 


tion mechanism. Further optical (0.5—0.6 fim) polarimetric 
studies (Capetti et al. 19951 using the Hubble Space Tele¬ 


scope (HST) showed a centrosymetric polarisation pattern 
along the ionisation cones in the ~10” (600 pc) central re¬ 
gion. This polarisation pattern is the signature of a cen¬ 
tral point source whose radiation (which can be polarised 
from multiple scattering into the funnel of the torus and/or 
broad line clouds) is ultimately polarised through scattering 
by dust and/or electrons. At ultraviolet and optical wave¬ 
lengths, these studies showed that electron scattering is the 
dominant polarisation mechanism in the core and ionisation 
cones of NGC1068. In the NIR, an increase in the polari¬ 
sation flux density with increasing wavelength in the core 


was observed (e.g. Lebofsky, Kemp fc Rieke 1978 Brindle 
|et al.|[r990| ). Further modeling using NIR polarimetric ob¬ 
servations has shown that the polarisation in the core most 
likely arises from dichroic absorption of radiation from the 
central engine by aligned dust grains in the torus ([Young et 


ar]|1995[ [Packham et al.||1997[ [Lumsden et al.||1999| |Sirnp 


son et al. 2002 Watanabe et al. 20031. A combination of 


electron and dust scattering from the ionisation cones still 
has some contribution, i.e. <10% of the polarised flux at 
2.0 pm ( Young et al.| 1995 1. Honig, Prieto fc Beckert (20081 
estimated that the contribution of the accretion disc emis¬ 
sion to the total flux is negligible at 2.0 pm, and that dust 
emission from a clumpy torus is the dominant contributor. 
In addition, [Watanabe et alT] ( |2003[ ) raise the possibility that 
scattering off large grains in the torus could also explain the 
NIR polarisation. They argued that the change of position 
angle (P.A.) of polarisation and the increase in the degree 
of polarisation from optical to NIR due to dust scattering 
by large grains in the torus cannot be ruled out. Although 
detailed studies have not been done yet, some implications 
are discussed in Section f4.21 

In the present paper, we examine the MHD outflow 
wind model in the core of NGC 1068 assuming the torus to 

||2008a|bp 


be clumpy (Nenkova et al. 2002 Nenkova et al. 


We follow the approach from Lopez-Rodriguez et al.|(|2013|) 
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The torus of NGC 1068 as a hydromagnetic outflow wind 3 


to estimate the magnetic field strength in the torus of NGC 
1068 and extent the stndy to estimate physical parameters 
of the torus using MHD outflow wind models, allowing us 
to describe the evolution and kinematics of the clumps. We 
performed imaging polarimetric observations in the J' and 
K' filters using MMT-Pol in conjunction with the adaptive 
optics (AO) system on the 6.5-m MMT. The paper is or¬ 
ganised as follows: Section describes the observations and 
data reduction, Section presents our polarimetric results, 
which are analyzed in Section]^ In Section[^we estimate the 
magnetic field strength in the torus of NGC 1068. Sectionj^ 
presents the discussion and Sectionthe conclusions. 


Table 1. Summary of Observations. 


Object 

Filter 

Frame Time 

(s) 

# ABA 

Total Time 

(s) 

NGC 1068 

J' 

15 

3 

540 


K' 

10 

10 

1200 

HD 38563C 

J' 

10 

2 

240 


K' 

10 

3 

360 


2 OBSERVATIONS AND DATA REDUCTION 

NGC 1068 was observed on 2013 October 23 using MMT- 
Pol ( [Packham fc Jones|2008[ [Packham et al.|2010[ |2012p in 
conjunction with the AO system and the f/15 camera on 
the 6.5-m MMT, Arizona. The AO system was used to as¬ 
sist the polarimetric observations of MMT-Pol. Specifically, 
a dichroic at 15° to the normal before the aperture window 
of MMT-Pol reflects optical light up to a CCD-based wave- 
front sensor and passes the IR beam into MMT-Pol at the 
Cassegrain focus. Thus, NGC 1068 is used for the AO cor¬ 
rection at optical wavelengths. MMT-Pol uses a 1024 x 1024 
pixels HgCdTe AR-coated Virgo array, with a pixel scale of 
0”043 pixel” corresponding to a field of view (FOV) of 
44 " X 44". MMT-Pol uses a rectangular focal plane aper¬ 
ture, a half-wave retarder (half wave plate, HWP), and one 
of two Wollaston prisms. The rectangular focal plane aper¬ 
ture provides two non-overlapping rectangular images with 
an individual FOV of 20" x 40". In standard polarimetric 
observations, the HWP is rotated to four P.A. in the follow¬ 
ing sequence: 0°, 45°, 22?5, and 67?5. A Calcite Wollaston 
is used in the wavelength range of 1—2 ^m and a Rutile 
Wollaston is used in the wavelength range of 2—5 /rm. 

The J' (Ac = 1.33 /rm, AA=0.08 /rm, 50% cut-on/off) 
and K' (Ac=2.20 ^m, AA=0.11 fj,m, 50% cut-on/off) filters 
provide the best combination of sensitivity and wavelength 
range within the instrumental filter set; thus, these filters 
were used for the observations. In both filters, the images 
were acquired in an ABA dither pattern with an offset of 
8" in declination, where four HWP P.A. were taken in each 
dither position. The position of the short axis of the array 
with the north on the sky was 149° and 46° E of N at J' 
and K', respectively. At J', frame exposure times of 15s per 
HWP P.A. at each dither position were taken, with a total 
of 3 ABA dither patterns, providing a total exposure time of 
540s. At K', frame exposure times of 10s per HWP at each 
dither position were taken, with a total of 10 ABA dither 
patterns, providing a total exposure of 1200s. A summary 
of observations is shown in Tabled] 

The data were reduced using custom IDL routines and 
standard NIR imaging procedures. The difference for each 
correlated double sample (CDS) pair was calculated, and 
then sky subtracted using the closest dither position in time 
to create a single image per HWP P.A. For each dither po¬ 
sition, the ordinary (o-ray) and extraordinary (e-ray) rays, 
produced by the Wollaston prism, were extracted and then 
the Stokes parameters, 7, Q, and U were estimated according 
to the ratio method (e.g. [Tinbergen 20061 for each dither po¬ 


sition. The Stokes parameters were registered and shifted to 
a common position, and then co-averaged to obtain the final 
7, Q, and U images. Finally, the degree, P = \/+ 77^/7, 
and P.A., PA = 0.5arctan(7//(3), of polarisation were es¬ 
timated. During this process, individual photometric and 
polarimetric measurements were performed for each dither 
position, allowing to examine for a high and/or variable 
background that could indicate the presence of clouds or 
electronic problems. Fortunately, no data needed to be re¬ 
moved for these reasons, but some data were removed when 
the AO guide unlocked. The measurements of the degree of 
polarisation were corrected for polarisation bias using the 
approach by Wardle &i Kronberg (19741. 


The image of the unpolarised standard star, HD 224467, 
in the K' filter was used to estimate the image quality of 
the observations, ensuring the quality of the AO system 
prior the science observations. For the standard star, the full 
width at half maximum (FWHM) using a Gaussian profile 
was estimated to be 0"17 x 0"23 in the K' filter (Figure]^ 
bottom-middle). MMT-Pol has a residual astigmatism that 
produces an X-shape on the FWHM. The polarised stan¬ 
dard star, HD 38563C, was observed in both filters. The ob¬ 
servations were performed during windy conditions on the 
observatory. The unpolarised standard star, HD 224467, was 
observed against the wind, minimizing instrumental shakes, 
whilst NGC 1068 observations were performed in favor of the 
wind, which affect the image quality. Frame exposure times 
of 10s per HWP P.A. at each dither position were taken, 
with a total of 2 and 3 ABA dither patterns, providing a 
total exposure time of 240s and 360s in the J' and K' filters, 
respectively. A summary of the observations is shown in Ta¬ 
ble The observations show HD 38563C as an elongation 
due to the nebulae that the star is embedded. Although the 
FWHM of the observations were affected, the polarimetric 
measurements are not affected as these measurements were 
estimated in a large circular aperture assuming HD 38563C 
as a single source. In polarimetry, HD 38563C allows us to 
estimate the zero-angle calibration of the observations. The 
zero-angle calibration, A9, was estimated as the difference 
of the measured P.A. of polarisation of our observations, Oji 
= 25°± 5° and 9^i = 37°± 2°, and the P.A. of polarisation, 
9Y = 71°± 1° and 9]^ = 78°± 17° provided by |Whittet ^ 


al. (19921. Thus, the zero-angle calibrations were estimated 
to be A9y = 9T - 9y = 46° ± 5° and A6 »k' = C - 9^, = 
41°± 17°. The instrumental polarisation was measured to 
be 0.05% ± 0.03% using several unpolarised standard stars 


from Whittet et al. (19921 during the observing run. 
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3 RESULTS 


Figure shows the J' (top-left) and K' (bottom-left) total 
flux images of NGC 1068. The extension towards the West 
in the K' total intensity image (Figure [^black diamond) 
is a ghost which we attribute to faint, internal reflections 
within the optics. This faint ghost does not appear in po¬ 
larised intensity. The radial profiles of the total flux images 
at J' (top-right) and K' (bottom-right) of NGC 1068 and 
the unpolarised standard star, HD224467, at K' are shown. 
The FWHMs in the longest and shortest direction and as¬ 
sociated position angle, using a 2D Gaussian profile, were 
estimated to be l'.'23 x I'.'Ol at P.A. = 56° E of N, and 
0'.'48 X 0'.'38 at P.A. = 50° E of N in the J' and K' fil¬ 
ters, respectively. Note that the FWHMs in Figure are 
calculated using the radial profiles. The estimated FWHM 
of NGC 1068 differs from those (~0Vl2 at 2.2 /rm) previ¬ 
ously reported in the literature using high-spatial resolution 


observations (e.g. Rouan et al. 1998 2004 Gratadour et 
al.||2006 ). Although the FWHM of NGC 1068 was affected 


during the observations, only those observations with (1) 
photometric conditions; (2) locked AO in a complete ABA; 
and (3) good image quality, were used. Specifically, as the 
nucleus of NGC 1068 is faint, the AO system was forced to 
sample the optical image to lower frequencies compared to a 
bright star. At these lower frequencies, the AO system com¬ 
pensate less for seeing, exacerbated by the wind buffeting. 
To account for the impact of the weather condition and low 
AO camera frame rate in the measured flux, we estimated 
the variation of the counts between the several ABA’s to be 
12%. This result is included in the uncertainties of the AGN 
flux contribution in Section [4.II 


Figure shows the polarised flux image and polarisa¬ 
tion vectors through the K' filter. The overlaid polarisation 
vectors are proportional in length to the degree of polar¬ 
isation with their orientation showing the P.A. of polari¬ 
sation. To improve the signal-to-noise ratio (S/N) and to 
obtain statistically independent polarisation measurements, 
the normalised Stokes parameters q = Q/1 and u — U/1 
were binned from 4x4 pixels (0'.'l72x0'.'l72) to a single 
pixel. Then, only those polarisation vectors with P/up > 3 
are shown in Figure Up is the uncertainty on the de¬ 
gree of polarisation within the binned pixels, estimated as 
= x/cTp -f cTp, where and are the uncertainties on 
the normalised Stokes parameters. In polarised flux, NGC 
1068 shows a resolved nucleus with a FWHM estimated to 
be 0'.'49 X 0'.'45 at P.A. = 30° E of N. The polarisation vec¬ 
tors and polarised flux image in the J' filter are not shown 
because of the low S/N ~ 40 at the peak pixel of the obser¬ 
vations. 


To compare with previously published values (Table [^, 
we made measurements of the nuclear polarisation at J' and 
K' in several apertures. In all cases, polarimetric errors were 
estimated by variation of the measurements in subsets of the 
data. Our measured polarisation values in a 2" aperture are 
in close agreement with previously published values (Table 
[^. [Simpson et J) (20021 measured a degree of polarisation 
of 6.0% ± 1.2%, with variations of ±2-3% in the degree of 
polarisation between the two observing runs of NGC 1068, 
however the reason of this difference is unknown. We sub¬ 
sequently conclude that the measured polarisation of |Simp-] 
son et al. (20021 and ours, 4.8% ± 0.1%, are in marginal 


Table 2. Nuclear Polarimetry of NGC 1068 and Literature Val¬ 
ues. 


Aperture 

(") 

Filter 

P 

(%) 

P.A. 

(°) 

Ref(s). 

0.2 

2.0 fim 

6.0 ± 1.2 

122 ± 15 

a 


K' 

4.8 ± 0.1 

127 ± 2 

b 

0.5 

J' 

2.9 ± 1.3 

111 ± 14 

b 


K' 

4.4 ± 0.1 

127 ± 2 

b 

2.0 

J 

2.25 ± 0.25 

106.8 ± 2.0 

c 


J 

1.71 ± 0.22 

106.45 ± 1.70 

d 


J' 

2.4 ± 1.2 

106 ± 17 

b 


K„ 

4.57 ± 0.50 

120.2 ± 2.0 

c 


K„ 

4.11 ± 0.46 

120.55 ± 2.38 

d 


K' 

4.2 ± 0.3 

119 ± 5 

b 


References; (a) [Simpson et aL]||2002[|; (b ) This work; (c) |Lumsden| 
jet al.j ( [l999| l; (d) jPackham et al.| ( |1997t 


agreement within the uncertainties. We found an increase 
in the degree of polarisation with increasing wavelength, as 
expected from previous polarimetric studies (e.g. Lebofsk^ 


Kemp fc Rieke|1978||Brindle et al.|1990||Young et al.|1995| 


Packham et al.| 1997 Lumsden et al.||1999 1 with a roughly 

constant P.A. of polarisation in the NIR wavelength. 

We also made measurements of the nuclear polarisation 
at J' and K' in a 0”5 (30 pc) aperture. These measurements 
minimize contamination from extended (diffuse) dust emis¬ 
sion, as well as obtain a statistically independent measure¬ 
ment with higher S/N than the single polarisation vectors 
(important in the J' filter due to the low S/N of the observa¬ 
tions). Specifically, its degree of polarisation was measured 
to be 2.9% ± 1.3% and 4.4% ± 0.1% with a P.A. of polarisa¬ 
tion of 111°± 14° and 127°± 2° in a 0V5 (30 pc) aperture in 
the J' and K' filters, respectively. The statistical significance 
of the degree of polarisation measurements is 2.4cr and 36cr 
in the J' and K' filters, respectively. 

The polarised flux image at K' shows a centrosym- 
metric polarisation pattern (SouthWest Lobe in Eigure 
at ~l'.'5 South from the center of the galaxy. The degree 
of polarisation is 15% ± 2% within the binned 4x4 pix¬ 
els (0” 172x0” 172). This polarisation pattern is the signa¬ 
ture of a central source whose radiation is scattered by 
dust/electrons. This polarised structure and measured po¬ 
larisation are in excellent agreement with (1) the 2.0 ^m po¬ 
larimetric observations using HST/NIGMOS by Simpson et 


al. (20021, and (2) the ionisation cones, with opening angle 


of 65°±20° and P.A.~15°, observed using HST [OHI] obser¬ 
vations by Evans et al.| ( [T9M| . We found a ~49° misalign¬ 
ment between the measured P.A. of polarisation at K' and 
the radio jet axis (Figure]^ insert plot). This misalignment 
is similar to the ~45° misalignment between the compact 
(1.35 X 0.45 pc) hot (~800 K) dust component (i.e. obscur¬ 


ing torus) observed by Raban et al. (20091 and the radio jet 


axis. A further interpretation of this misalignment in terms 
of the MHD outflow wind model in AGN is given in Section 

IQ 
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NGC 1068 
FWHM = 1.12" 



Offset (") 


0,6 0,8 1.0 
Offset (") 




Offset (") 


Offset 


Offset (") 


Figure 1. The 8”5 x 8”5 central region of NGC 1068 in the J' (first row) and K' (second row) total flux (grey scale) images. In both 
images, the lowest level contour is at So-, then contours are at the sigma-levels of 2", where n = 4, 5, 6, etc. The J' image was rotated 
for direct comparison with K' image (Section]^. The ghost (black diamond) in the K' is marked. The middle plot of the second row 
shows the unpolarised standard star, HD224467, at K' with the same FOV and contour levels as the total flux images of NGC 1068. In 
all the total flux images, the pixel scale is 0.043” pixel“^. The radial profiles of NGC 1068 (black solid line) and the standard star (red 
dashed line), normalizes to the flux density value of the centroid of NGC 1068, are shown. 


4 ANALYSIS 

4.1 AGN contribution to the total flux density 

The emission from the unresolved core of NGC 1068 at K' 
can arise from two components: (1) diffuse stellar emission in 
the nuclear bulge, and (2) AGN emission. We interpret the 
AGN emission as a point-source emission component, while 
the diffuse stellar emission is extended. Photometric cuts 
through the nucleus in the J' filter show negligible evidence 
of a nuclear point-source, AGN emission, due to extinction 
towards the core and/or contamination by the diffuse stellar 
emission (Figure]^ top-right). The J' filter profile is assumed 
to be representative of the diffuse stellar emission in the nu¬ 
clear bulge. At K', both the AGN emission and diffuse stellar 
emission are detectable. To estimate the relative contribu¬ 
tions from both emission components in the K' filter, two 
different methodologies were followed. 

In the first method (method 1), we followed similar anal- 


to that of Turner et al. 

(1992 

; Packham et al. (19961 

Lopez-Rodriguez et al. 

(2013 

I. We took photometric 


profiles along the major axis of the diffuse extended stellar 
emission in the J' and K' filters, which is shown to have a 
constant extinction of Ay = 0.7 ± 0.7 ( Young et al.|[lOT5 |. 
These cuts allow us to characterize the host galaxy contribu¬ 
tion. The contribution of the AGN emission is assumed to be 
the point-spread function (PSF) of the observations (i.e. the 


unpolarised standard star, HD224467). To fit the K'-profile, 
the emission was modeled as the sum of a scaled-PSF and 
the diffuse stellar emission from the J'-profile. Using this 
method the best estimate of AGN contribution to the total 
flux is 60% ± 12% in a 075 (30 pc) aperture. The uncer¬ 
tainty was estimated as the 12% uncertainty due to weather 
conditions and low AO camera frame rate (Section]^. 

In the second method (method 2), the PSF from method 
1 was scaled to the flux density of the K' emission to subtract 
the AGN emission. This procedure produces a “flat-top” 
profile over the unresolved core of NGG 1068 (for examples 
of the use of this technique, see|Radomski et al. |2002 |2003| 


Ramos Almeida et al. 2011 Lopez-Rodriguez et al. 


20131. 


Using this method, the contribution of the AGN emission to 
the total flux is 50% ± 19% in a 075 (30 pc) aperture. The 
uncertainty was estimated as the addition in quadrature of 
the uncertainty of the “flat-top” method applied, 15%, and 
the 12% uncertainty due to weather conditions and low AO 
camera frame rate (Section]^. 

From both methods presented above, the average of the 
AGN emission to the total flux is estimated to be 55% ± 
16%, whereas the contribution of the diffuse stellar emission 
to the total flux is 45% ± 16%, in a 075 (30 pc) aperture in 
the K' filter. The formal total uncertainties are 16% through 
the average of both methods. Systemic errors of methodolo¬ 
gies could increase these estimates, but are difficult to quan¬ 


tify. For NGG 1068, Weinberger et al. (19991 estimated the 
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Offsets (") 

Figure 2. The 8" X 4" polarised flux image (grey scale) at K' with overlaid polarisation vectors. The contour levels are 1.5%, 5%, 
10%, 30%, 50%, 70% and 90% of the peak of the polarised flux image. The polarised flux image has a pixel scale of 0”043, whilst the 
polarisation vectors with PfcTp > 3 within the binned 4x4 pixels box (0^172) are shown. A vector of 5% of polarisation and the FWHM 
of the observations are shown. The dashed line shows the orientation of the radio jet at P.A.~ 12° by|Gallimore et al.| | |1996[ l, and the 
yellow shadow area shows the HST [OIII] opening angle, 65°ib20° with P.A.~15°, by [Evans et al.| | |1991| . The P.A. of polarisation and 
the radio jet axes (insert plot) with respect to the North is shown (Section]^. 


AGN contribution to be 49% ± 8% of the total nuclear flux 
in a 0('48 aperture in the K filter using the near infrared 
camera (NIRC) on the 10-m Keck I telescope. jWeigelt et al.j 
(20041 estimated a total flux contribution of 46% ± 5% at 
K' in a 18 X 39 mas using NIR speckle interferometry on 
the 6-m Special Astrophysics Observatory (SAO). 


4.2 The nuclear intrinsic polarisation 

The diffuse stellar emission in the nuclear regions of NGC 
1068 significantly dilutes the observed polarisation from the 
AGN. If only the dilution from the diffuse stellar emission 
is subtracted from the observed polarisation, the estimated 
intrinsic polarisation is independent of the dominant po¬ 
larisation mechanism in the nuclear regions of NGC 1068. 
Thus, the dominant polarisation mechanism in the K' filter 
is needed to estimate the intrinsic degree of polarisation of 
the AGN. The intrinsic degree of polarisation arising from 
the dominant mechanism of polarisation in the nuclear re¬ 
gions of NGC 1068 is estimated as follows. The measured 
polarisation at K', = 4.4% ± 0.1%, in a 0”5 (30 pc) 

aperture is corrected by accounting for (1) the measured de¬ 
gree of polarisation through dichroic absorption of starlight 
in an off-nuclear region (P^/^), (2) the diffuse stellar emis¬ 
sion relative to the AGN emission in the nuclear aperture 
and (3) the dominant polarisation mechanism to the 
polarised flux of the AGN in the nuclear aperture (P^/°). 

The degree of polarisation at K' through dichroic ab¬ 
sorption of starlight in the off-nuclear regions of NGC 1068 
was measured to be P^*7 = 0.5% ± 0.3%. The off-nuclear 


polarisation was estimated as the average of several polari¬ 
sation measurements using a 2” aperture. We assume that 
the off-nuclear polarisation has the same P.A. of polarisation 
as the AGN polarisation. Any value of the P.A. of polarisa¬ 
tion will produce negligible variations in the estimation of 
the intrinsic polarisation, given the low measured polarisa¬ 
tion. Using the diffuse stellar emission, = 45% ± 16%, 
relative to the AGN emission, = 55% ± 16% estimated 
in Section |4.1[ we define the ratio of both emission compo¬ 
nents as Pk' = = 0.8±0.6. 


As noted in the introduction, previous studies (Young et 


aI.|1995||Packham et al.|1997||Lumsden et al.|1999||Simps~ 


et al.|2002 Watanabe et al.|2C103 \ have shown that the NIR 

polarisation in the nuclear regions of NGC 1068 arises from 
the passage of light through aligned dust grains in the torus. 
As we only have polarisation measurements in two NIR fil¬ 
ters, we cannot estimate the relative contribution of the 
dichroic absorption to the polarised flux at K'. We therefore 


resort to previously published results. Young et al. (19951 


found that the dichroic absorption mechanism accounts for 
~90% of the polarised flux, with a combination of electron 
and dust scattering from the ionisation cones accounting for 
~10% of the polarised flux at 2.0 /rm in a 3'.'08 x 3700 (185 
X 180 pc) slit aperture (figure 6 from Young et al. 19951. 
Watanabe et al. ( 2003| ) found that the dichroic absorption 
mechanism can account for ~100% of the polarised flux at 
2.2 /rm in a 0788 x 2700 (53 x 120 pc) slit aperture (fig¬ 
ure 1 there). These slit sizes are larger than the 0('5 (30 pc) 
aperture in our study. The difference in aperture sizes can 
affect the polarisation measurements due to scattered emis¬ 
sion from the extended biconical structures around the core 
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The torus of NGC 1068 as a hydromagnetic outflow wind 7 


of NGC 1068 (Packham et al. 1997 Lumsden et al. 1999 
Simpson et al.|2002 |, and the diffuse stellar emission in the 
host galaxy. [Simpson et ^ ( |2002[ ) showed that the polarised 
flux from the core of NGC 1068 completely dominates within 
an aperture <l", where contributions from off-nuclear scat¬ 
tered light would be barely detectable. Other polarisation 
components (i.e. electron scattering from the inner regions 
of the ionisation cones and/or dust scattering by large grains 
in the torus) can contribute, at some level, to the polarised 
flux at K' within our 0”5 (30 pc) aperture, however the con¬ 
tribution of these components are highly dependent of the 
geometry of the system, implying a large amount of assump¬ 
tions that will introduce high uncertainties in the estimation 
of the intrinsic polarisation. Based on these studies, these 
polarisation mechanisms cannot be ruled out, but the con¬ 
tribution to the polarised flux of these components have not 
been quantified yet and require more detailed polarisation 
models. We here assume the upper limit of the dichroic ab¬ 
sorption to be 100 % of the polarised flux, in agreement with 
the polarisation model of Watanabe et al. (20031, i.e. 

= 1. The intrinsic polarisation at K' arising from dichroic 
absorption, Pfl, is 


P'fl + 


( 1 ) 


and is estimated to be P/T = 7.0% ± 2.2% in a 0V5 (30 pc) 
aperture. 


5 MAGNETIC FIELD STRENGTH IN THE 
TORUS 

The magnetic field strength in the NIR emitting regions of 
the torus of NGG 1068 is estimated following similar analy- 


K. The dust grain sizes are assumed to be in the range of 
0.005—0.25 /rm. Based on the Clumpy torus description and 
the further fitting to the nuclear IR SED of NGC 1068 by 


Alonso-Herrero et al. (20111 using BayesClumpy (Asensio 


Ramos & Ramos Almeida 20091, the torus of NGC 1068 
can be described with a radial thickness of E = 6 /li, angu¬ 
lar width of CT = 26 ^ 4 °, number of clouds along the equa¬ 
torial No = Idlg, radial density profile with an index of 
q = 2 . 21 q 3 , viewing angle of i = 88 ^ 3 ° (almost edge-on), 
and the optical depth per single cloud of tv = dOlg. The 
distance of the inner wall of the torus from the black hole, 



and we estimated rs^b = 0.4 pc for dust grains at 1500 K. 
The torus height is estimated as 77 = rout sin cr = l.llg'g 
pc, with Tout = Vr^ub = 2-3.2 pc. The half-opening angle 
of the torus is estimated as 0 = 90° — a = 64lg°, which 
is defined with respect to the projected symmetry axis of 
the Clumpy model. The number density in the individual 
clumps of NGC 1068 were calculated assuming the output 
parameters of the Clumpy torus models by | Alonso-Herrei^ 
|et al.| ( |2011[ ): (1) a torus radius of rout = 2.0-3.2 pc, (2) 
the number of clouds along the equatorial direction of No = 
I 4 I 3 , and (3) the optical depth per cloud of Tv = dOlg, con¬ 
verted to be Av = 1.086rv = 53 I 3 mag. Using the s tandard 
Galactic ratio A^/Nh = 5.23 x 


10 


mag cm (|Bohlin, 


Savage & Drakell978 1 , the column density of the individual 

_ . I no . _ no n_ 


clouds is Ah = i.Ol/jgQgX 10^^ cm“^. The number density 
in the individual clouds in the torus are estimated to be n = 
Nn/iraut/No) = 2.30^0.61 ^ 10 ® cm“®, which compare very 
well with the number density in the range of 10 ^ — 10 ® cm“® 
for molecular clouds in Orion A, M17 and Cepheus A (e.g. 


sis to that of |Lopez- Rodriguez et al.|(|201^. Specifically, the |Goldsmith[|1999[). A summary of the physical parameters of 


magnetic field strength is estimated through three different 
methods: ( 1 ) paramagnetic alignment, ( 2 ) thermal and mag¬ 
netic relaxation time equipartition, and (3) Chandrasekhar- 
Fermi method. 


the torus is shown in Tabled 


5.1 Physical conditions in the torus of NGC 1068 

NIR reverberation mapping of several AGNs has shown that 
the outer radius of the broad line region (BLR) approxi¬ 
mately corresponds to the inner radius of the dusty torus 
( Suganuma et aL]|2006 Koshida et al.|[2014 |. The gas tem¬ 
perature reaches a value of ^10“* K in the BLR (Netzer 


1987). Krolik & Kriss (20011 suggested that a warm ab¬ 


sorber gas in the inner edge of the torus can reach temper¬ 
atures in the range of 10^ — 10® K. Based on these previous 
studies, we adopt the lower-limit of the gas temperature to 
be Tgas = 10"^ K. We here adopt the Clumpy torus modeQ 
by Nenkova et al. (2002); Nenkova et al. ( 2008a|b ). This 
model depicts the torus as a distribution of optically thick 
and dusty clouds surrounding the central engine, instead of 
homogeneously filling the torus volume. The NIR-emitting 
region in the torus is identified with clumps directly illumi¬ 
nated by the central engine, with dust grain temperatures 
below the sublimation temperature in the range of 800—1500 


^ For details of the Clumpy torus models see www.clumpy.org 


5.2 Method 1: Paramagnetic alignment 

We concluded that the NIR polarisation arises from the pas¬ 
sage of radiation through aligned dust grains in the clumps 
of the torus of NGC 1068 (Section |4.2[ ). Non-spherical spin¬ 
ning dust grains can become aligned in the presence of a 
magnetic field, a process termed paramagnetic alignment 
( Davis &; Greenstein||195T l. If the dust grains are aligned by 
this mechanism, then some degree of polarisation through 
dichroic absorption and/or emission can be measured. This 
measurement can be used to estimate the strength and ge¬ 
ometry of the magnetic field. In case of the geometry, the 
long axes of the dust grains are aligned perpendicular to 
the direction of the magnetic field. A P.A. of polarisation 
parallel to the magnetic field is expected for dichroic ab¬ 
sorption. In case of the magnetic field strength, if paramag¬ 
netic alignment is the dominant mechanism, then the dust 
grain alignment efficiency is given by the ratio of the de¬ 
gree of polarisation, P(%), and the extinction, Av, and is 
related to the magnetic field strength, B, as sho wn by | Jones | 
& Spitzer (1967), and adapted by Vrba et al. ( 1981L and 


Lopez-Rodriguez et al. (2013): 
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8 E. Lopez-Rodriguez et al. 


Table 3. Physical Parameters of the Torus of NGC 1068 


Description 

Parameter 


Value 

Radial extent of the torus 

Y 


ell 

Width of the angular distribution of clouds 

<J 


26l^° 

Number of clouds along the radial equatorial direction 

No 


1411 

Power-law index of the radial density profile 

Q 


0 q+0-4 

^•^-0.3 

Inclination of the torus 

i 


88l3° 

Optical depth per single cloud 

Tv 


49l^ 

Sublimation radius 

Tsub 


0.4 pc 

Torus height 

H 


I+I 0.3 P= 

Half-opening angle 

e 


64l^° 

Gas temperature 

T 

J- gas 


104 K 

Grain temperature 

T 

Tgr 

800 

- 1500 K 

Grain size 

a 

0.005 

— 0.25 pm 

Number density in the clumps 

n 

2.30l°;|® X 

10 ® cm”® 




1/2 


(7- 1) 


T 


1 - ^ ) ( 2 ) 


x" is the imaginary part of the complex electric susceptibil¬ 
ity, a measure of the attenuation of the wave caused by both 
absorption and scattering, B is the magnetic field strength, 
7 is the ratio of inertia momentum around the two major 
axes of the dust grains, Tg, is the dust grain temperature, a 
is the dust grain size, w is the orbital frequency of the dust 
grains, n is the number density in the cloud, mu is the mass 
of a hydrogen atom, k is the Boltzmann constant, and Tg^a 
is the gas temperature. 

Davis & Greenstein ( 1951| showed that the lower bound 


of the ratio x for most interstellar grains is: 


^ = 2.5 X 10"^^ T~R 


(3) 


The ratio of the moments of inertia of the dust grains, 
7 , is defined as: 


1 

^=2 


- 1+1 
a 


(4) 


where 6/a is the dust grain axial ratio. A typical value of 
6 /a for interstellar dust grains is ~ 0.2 (Aannestad & Purcell 
1973[ Kim & Martin 19951. It should be noted that the real 


dust composition in AGN might be more complex that the 


typical interstellar dust g 

;rains, as several authors suggest 

(e.g. Maiolino et al. |2001 

Maiolino, Marconi & Oliva|2001 

Hatziminaoglou E. et al. 

20151 


The magnetic field strength is estimated using the phys¬ 
ical conditions in Table the intrinsic polarisation at K' 
arising from dichroic absorption of P//* = 7.0% ± 2.2% (Sec¬ 
tion |T2|^_and_the_extinction to the central engine Av = 36 
mag (Packham et al.||l997l convertec0to A^' = 4.03 mag. 


The extinction represents the closest measurement at our 
observations, and it has been previously used by |Watanabe| 


The conversion of visual extinction. Ay, to extinction at K, Ak 


is Ak = 0.112Av I Jones|1989^. 


|et al.| ( [^003[ ). |Imanishi et al.| ( |1997[ ) estimated a visual extinc¬ 
tion in the range of 17—30 mag using the 3.4 /rm absorption 
feature of the NIR (2.9—4.1 /rm) spectroscopic observations 
on the 1.5-m telescope on Mount Lemmon, Arizona, USA. 
These measured visual extinctions are lower than the whole 
visual extinction of the clouds, Av = 53 I 3 mag. (Section 
EU, estimated using the Clumpy model to the IR SED of 
NGC 1068, indicating that we are possibly seeing the radia¬ 
tion from the clouds and/or central engine passing through 
a section of the clouds (Section |6.1[ ). Based on the range of 
the physical parameters shown in Table^ the magnetic field 
strength is estimated to be in the range of 4—82 mG for the 
NIR emitting regions in the torus of NGC 1068. The mag¬ 
netic field strength in the maser emission of NGC 1068 was 
estimated to be 71 mG and 3 mG for standard accretion disc 
and advection-dominated accretion flow (ADAF) model, re- 


spectively | 

Gnedin et al. 

2014 

1 . Several studies (e.g. 

Her- 

rnstein et al.||19981 |Kartje 

, Konigl & Elitzurll 19991 IModjaz 

et al. 20051 have shown that the magnetic field strength 


is in the range of 2—300 mG at scales of 0.1—1 pc from 
the central engine of AGN. Although these measurements 
were done for different objects, our estimated magnetic field 
strength in the torus of NGC 1068 is in agreement with the 
typical magnetic field strengths in the surrounding areas of 
the accretion disc in AGNs. 


Several assumptions were made in the estimation of the 
magnetic field strength in the torus of NGC 1068. Here, 
we consider each of these assumptions in detail. The phys¬ 
ical conditions in molecular clouds make the alignment of 


dust grains less responsive to the magnetic field (i.e. Lazar- 


1995 Gerakines et al. 19951. However, the dust grains 


in the dusty torus of AGN certainly experience much more 
turbulence and extreme physical conditions than in molec¬ 
ular clouds, and the dust alignment is more sensitive to the 
magnetic field strength (Section 5.3 6.11. 


We assumed a homogeneous magnetic field, where any 
inhomogeneities of the magnetic field in the torus are ig¬ 
nored. This assumption has some implications for the ratio 
P(%)/Av. If a homogeneous magnetic field is responsible 
for dust grains’ alignment in the torus, then all the dust 
grains will be aligned along the same direction of the mag¬ 
netic field. In this case, the alignment efficiency, P(%)/Av, 
is maximal. The degree of polarisation would decrease when 
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The torus of NGC 1068 as a hydromagnetic outflow wind 9 


inhomogeneities of the magnetic field are present (Section 


6.1 1 . This method shows a strong dependence on dust grain 


sizes. The grain size distribution in the clumps is poorly 
constrained, which makes it difficult to quantify the uncer¬ 
tainties introduced in the estimations of the magnetic field 
strength. 


5.3 Method 2: Thermal and magnetic relaxation 
time equipartition 


The approach followed in Method 1 (Section |5.2[ ) is valid 
only when the ratio of thermal to magnetic pressure is 
smaller than unity ( Gerakines et al.||1995 1, meaning that the 
magnetic field is dominant. Under this condition, the mag¬ 
netic field is strong enough to align the dust grains faster 
than the rotational kinetic energy. The magnetic relaxation 
time is required to be shorter than the thermal relaxation 
time, i.e. < Uti- To verify this condition, we calculate the 
required lower limit magnetic field strength given by (section 
7.2 from Lopez-Rodriguez et al.|[2013l: 


> 2.4 X lO^anmaTg, (' ( 5 ) 

V TrrriH ) 

The magnetic field strength is estimated to be >6 mG 
for the physical conditions shown in Table Under the con¬ 
dition that the magnetic relaxation time is shorter than the 
thermal relaxation time, the gas and dust are decoupled. 
This result has two implications: (1) it is only valid in low- 
density regions of the clouds, and (2) the ratio P{%)/A^ 
is dependent on the magnetic field strength. Thus, we con¬ 
clude that the detected NIR polarisation arises from the 
passage of NIR radiation through the low-density regions of 
the clumps in the torus. Otherwise, the NIR radiation would 
be completely extinguished (Section |6.1[ ). This scheme sat¬ 
isfies condition ( 1 ) and the estimation of the magnetic field 
strength in method 1 (Section |5.2[) is allowed. 


5.4 Method 3: Chandrasekhar-Fermi method 


The Chandrasekhar-Fermi method (Chandrasekhar & Fermi 


1953 hereafter CF method) gives us an alternative esti¬ 
mation of the magnetic field strength. Specifically, the CF 
method relates the magnetic field strength with the disper¬ 
sion in polarisation angles, a, of the constant component of 
the magnetic field, and the velocity dispersion of the gas. 
We use the adapted version of |Ostriker, Stone &: Gammiej 
(20011, where they introduced the factor of 0.5 in the CF 


method to compensate for averaging along the LOS. 


'A \ 1/2 

5 = 0.5 ( -^p\ ^ [pG] 


( 6 ) 


where p is the volume mass density in g cm~^, fjv is the 
velocity dispersion in cm s“^, and a is the dispersion of po¬ 
larisation angles in radians. The volume mass density was 
calculated using the number density in Table [3| mul tiplied 
by the weight of molecular hydrogen. Greenhill et al. (1996 1 
estimated velocity dispersions up to 100 km s“^ of the VLBI 
H 2 O masers at scales of 0.65 pc from the central engine of 
NGC 1068 using VLBI. We note that these results are higher 



Figure 3. Degree of polarisation at K versus optical depth. Plot¬ 
ted is the data compilation from pones| ( |198^ (black dots) with 
some specific objects labeled. The upper (long dash) line is the 
JKD model result assuming only a constant component to the 
magnetic field (P = tanh(rp)). The bottom (short dash) line 
corresponds to a purely random magnetic field with a decorrela¬ 
tion length of r/c = 0.1 (P oc r^/^). The solid blue line results 
from a 50/50 mix of constant and random components and fits 
the overall trend in the data well. The data point for IC 5063 
is plotted as a brown star | |Lopez-Rodrigue 2 et al^|2013[ |. The 
estimated intrinsic polaris ation for NGC 1068 fro m this work, 
P^J “ ^ 2.2% (Section 4.2 I with Ay = 36 from (Packham et 

|al.|l997| l, is plotted as a green star. 


than the typical velocity dispersion of 10 km s ^ used in pre¬ 


vious MHD outflow wind models (e.g. Elitzur & Shlosman 


20061 , which can yield higher estimations of the magnetic 


held strength. 

As the AGN in NGG 1068 is unresolved, the polarisa¬ 
tion vectors are average values of the LOS to the nucleus. If 
we average the LOS across the observed core (Figure]^, it 
will provide an insignihcant dispersion in P.A. This measure¬ 
ment will yield unrealistically large magnetic held strengths 
with the GF method. Hence, the dispersion of P.A., a, can¬ 
not be measured directly from the observed dispersion of the 
polarisation vectors across the core of NGC 1068 at K'. We 
need another approach to estimate a; we used the model by 


Jones, Klebe & Dickey (1992 hereafter JKD). This model 


relates the degree of polarisation at K with the level of turbu¬ 
lence in the interstellar medium and the magnetic held. The 
magnetic held is assumed to have a constant and a random 
component. For NGC 1068, using the intrinsic polarisation 
arising from dichroic absorption at K', 7.0% ± 2.2% 

(Section 4.2 (, and the extincti on to the central engine, = 
36 mag ( Packham et al.||l997 l, convertec^to be tk = 3.24, 
we found that our data in the JKD model is close to equal 
contribution of the constant and random components of the 
magnetic held (Figure]^. For our data point in Figure]^ 
if we assume that the constant component of the magnetic 
held is in the plane of the sky, the dispersion of polarisation 
angles is estimated to be a = 7.0° (0.1221 radians) using 


^ The conversion of visual extinction, Ak, 
Tk = 0.09Av I Jones|1989 I 


to optical depth at K, 
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Table 4. Estimations of the Magnetic Field Strength. 


Method 


B (mG) 


1; Paramagnetic alignment 4—82 

2: Thermal and magnetic equipartition >6 

3: Chandrasekhar-Fermi method 52 

4: Modified Chandrasekhar-Fermi method® 139 


+ 4: 


+ 11 
-20 


^Modified CF Method is calculated as B' = B X p with 
P = 0.15 (Section|6.1|l. 


figure 9 by Jones, Klebe & Dickey ( 1992| ). We substituted 
the above numerical values in Equation (1^, and estimated a 
lower limit of the magnetic field strength in the plane of the 
sky in the range of 52lg mG, depending on the conditions 
in the torus. 

We assumed that the constant component of the mag¬ 
netic field strength is in the plane of the sky. If the constant 
component of the magnetic field is away from the plane of 
the sky, then the magnetic field strength will be underesti¬ 
mated. For example, if the magnetic field is pointing along 
our LOS, zero polarisation will be measured. Since NGC1068 
is classified as a Type 2 AGN, the torus axis is in the plane 
of the sky. This geometry favors the assumption of the con¬ 
stant component of the magnetic field to be in the plane of 
the sky. We note that the all-sky Planck polarisation data 
( Planck Collaboration et al.|2014 l shows some objects with 
polarisation measurements above the theoretical polarisa¬ 
tion maximum in the JKD model. However, it is difficult to 
compare the Planck results with ours, given that Planck’s 
beam of 5' does not resolve AGNs. 

A summary of the estimations of the magnetic field 
strengths through the different methods presented in Section 
is shown in Table [4] The modified CF method described 
in Section ffi.II is also included. 


sured P.A. of polarisation and the radio jet axis (Section]^ 
Figure can be interpreted as the misalignment between 
the torus (assuming that the torus axis is in the plane of the 
sky) and the radio jet axis. Thus, physical components far 
below the spatial resolution of the NIR polarimetric obser¬ 
vations can be estimated. 

As noted in Section [5. 4[ we assumed the magnetic field 
as a mix of a turbulent and a constant component, where 
the turbulent component can be in any orientation. To infer 
the relative contribution of both components, we used the 
method presented by Hildebrand et al. |( |m^ . These authors 
showed that the ratio of turbulent to constant components 
of the magnetic field strength, /Bo, is a function of 

the dispersion of the plane of polarisation, a, given by: 


(B?) 


2\l/2 


Bo 




( 7 ) 


Using the dispersion of the plane of polarisation a = 


7° (0.1221 radians) from Section 5.4 the ratio of turbulent 


to constant components of the magnetic field strength is 
= 0.09. Thus, the constant component, with a 
toroidal geometry, is dominant in the torus of NGC 1068. 

The misalignments of the structures in the inner few 
parsecs of NGC 1068 can produce local inhomogeneities in 
the magnetic field. These inhomogeneities produce an under¬ 
estimation of the magnetic field strength measured through 


the paramagnetic alignment (Section 5.21. This underesti¬ 
mation of the magnetic field strength can explain the sim¬ 
ilarity of the paramagnetic alignment and the thermal and 
magnetic relaxation time equipartition. Thus, a higher mag¬ 
netic field strength should be expected. The CF method as¬ 
sumed thermal and magnetic pressure equipartition. |Kudoh| 
fc Basu[ (20031 showed that the equipartition is not satisfied 
in all regions of the molecular clouds and the CF method 
must be refined/modified in such regions. They defined the 
thermal to magnetic pressure ratio, /3, as: 


6 DISCUSSION 

6.1 Magnetic field in the torus of NGC 1068 


t-’ 


As discussed in Section [4.2| the measured polarisation arises 
from the passage of radiation through aligned dust grains in 
the torus of NGC 1068. As the grain alignment occurs with 
the long axes perpendicular to the magnetic field, a mea¬ 
sured P.A. of polarisation parallel to the magnetic field is 
expected for dichroic absorption. We interpret the measured 
P.A. of polarisation at K' in a 0'.'5 (30 pc) aperture, PA^' 
= 127°± 2°, as the orientation of the magnetic field in the 
torus of NGC 1068. We took the orientation of the torus to 
be PAtorua ~ —42° (measured anti-clockwise, for direct com¬ 
parison with our PAk', then PAtoma ~ —42° -|- 180°= 138°) 
from Raban et al. ( 2009| ), and we found that the PAtoma ~ 
PAp,'. The constant component of the magnetic field is thus 
aligned with the torus axis. Our result implies that our mea¬ 
sured P.A. of polarisation gives us direct information on the 
existence and geometry of the torus. Specifically, the polar¬ 
isation signature of the parsec-scale structure, identified as 
the obscuring torus, is observed through NIR polarisation 
and direct spatial comparison between IR interferometric 
and IR polarimetric observations can be made. Based on 
these results, the misalignment of ~49° between the mea- 


Ui (HV4^) 


( 8 ) 


where Cb is the sound speed defined as Ca = \JkTjm, k is 
the Boltzmann constant, T is the gas temperature, and m is 
the mean molecular mass. Va is the Alfven speed, defined as 
Va = Bjy/iirp, p is the cloud density, and B is the magnetic 
field strength. 

The dimensionless /3 parameter is typically used in 
MHD simulations (i.e. |Ostriker, Stone &i Gammie|2001[|Ku-| 
doh & Basu 20031 to account for the magnetization of the 


molecular clouds. Three different situations can be found: 
(1) (3 ^ 1 is considered as low magnetization, (2) /3 ~ 1 
represents the equipartition between dynamic and magnetic 
pressure (GF method condition), and (3) /3 ~ 0 is consid¬ 
ered as high magnetization. |Kudoh fc Basu| ( |2003[ ) found 
that the dynamic and magnetic pressure ratio is in approx¬ 
imate equipartition, /3 ~ 1, in the regions containing most 
of the mass of the molecular clouds. However, /3 < 1 in the 
low-density outer regions of the clouds. 

Using the magnetic field strength, B = 52lg mG from 
the GF method, and the cloud density from Section [5.4[ we 
find P ~ 0.15. This value suggest that the clumps in the 
torus of NGC 1068 are in a highly magnetic environment. 
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Figure 4. Sketch of a single cloud and its location within the 
torus (box). (1) The NIR total flux comes from the directly il¬ 
luminated clumps located above the equational plane affording a 
direct view into our LOS, i.e., the surface of the clumps that are 
directly illuminated by the central engine ( |Honig, Prieto fc Beck-| 
|ert||2008| l, (2) the NIR radiation passes through the low-density 
outer layers of the clumps, and (3) NIR polarisation arises from 
the passage of radiation through the aligned dust grains in the 
clumps (Section|6.1|l. 


Another implication of /3 ~ 0.15 is that the polarisation 
is arising from the low-density regions of the clumps. This 
physical condition agrees with the conditions of the param¬ 
agnetic alignment mechanism (Section |5.3[ ). Thus, the mag¬ 
netic field strength estimated through the CF method should 
be refined using the ft paramet er. The modified CF method 
is giv en hy B' = B x (Ostriker, Stone & Gammie 

20011, which yields a magnetic field strength B' = 139)!)2o 
mG (Table 1^. 

We interpret these results as: (1) the NIR total flrrx 
comes from the directly illuminated clumps located above 
the equatorial plane affording a direct view into our LOS, 
i.e., the surface of the clumps that are directly illuminated by 
the central engine (e.g. [Honig, Prieto fc Beckert|2008 Hey- 


mann fc Siebenmorgen||2012 Stalevski et al.||20i2 |, (2) the 

NIR radiation passes through the low-density outer layers of 
the clumps, and (3) NIR polarisation arises from the passage 
of radiation through the aligned dust grains in the clumps. 
These dust grains are aligned by a global toroidal magnetic 
field in the torus generated by the accretion disc. The clumps 
located above the equatorial plane, those clumps in (1), have 
higher likelihood of escaping photons along our LOS than 
those located close to the equatorial plane. A sketch is shown 
in Figure]^ 


6.2 Hydromagnetic wind 


Within the MHD framework (e.g. Blandford & Payne 1982 


Emmering, Blandford fc Shlosman|1992| |Gontopoulos 1994 


Kartje, Konigl fc Elitzur||1999 1, the magnetic field plays an 

important role in the dynamics of the clumps in the torus. 
The inflow/outflow mass rates for those clumps showing NIR 
dichroic absorption are estimated assuming a hydromagnetic 
wind model. The mass accretion rate as a function of the 
magnetic field strength can be estimated through the ther¬ 
mal and magnetic pressure ratio. The thermal pressure on 
the midplane (Frank, King fc Raine|[T9^l is given by: 


2 GMbuM 
Pth = pc„ = 


‘S'K'/Tkc 


( 9 ) 


where Mbh is the black hole mass, M is the mass accretion 
rate in units of Mq yr“^, r is the distance from the black hole 
in units of parsecs, and is the Shakura-Sunyaev ( [Shakur^ 
fc Sunyaev]1973 1 viscosity parameter. 


As described in Section [6.1[ the thermal and magnetic 
pressure equipartition was defined as ft = Pth/P- b, which can 
be used in Equation ||^ to estimate the mass accretion rate 
as a function of the magnetic field: 


M = 


3-s/^Cs r ® Qs/3 

SGMbh 


( 10 ) 


The parameters Os and ft are related through the MHD 
simulations prescription = G/{I -\- ft), where C is a con¬ 


stant typically set in the range of 0.5—0.6 (Narayan et al. 


19981. This relation is satisfied if magnetic fields provide 


the kinematic viscosity, which is the case in the magnetic 
environment in the clumps, ft ~ 0.15, estimated in Section 
|6.1| In this context, the viscosity parameter is estimated to 
be = 0.48 ± 0.05. We use the magnetic field strength 
estimated to be in the range of B = 4—139 mG (Section 
[^. We took the black hole mass of NGG 1068 Mbh = 
(8.0 ± 0.3) X 10® Mq from Lodato fc Berlin (20031. These 
authors fitted a self-gravitating accretion disc model to the 
observed non-Keplerian rotational curve of the maser disc to 
derive the black hole mass in NGC 1068. From the physical 
interpretation presented in Section [6.11 the NIR emitting re¬ 
gion arises from the directly illuminated faces of the clumps 
located as = 0.4 pc for dust grains at 1500 K (Section 


5.11. 


The clumps with measured NIR dichroism were as¬ 
sumed to be above the equatorial plane (Figure [^, these 
clumps are located in the inner edge of the torus, with a 
vertical heigh of B = r^ub sin cr = 0.17)(o'o2 PC. As a first 
approximation, the balance of magnetic and thermal pres¬ 
sure on the midplane allow us to estimate upper limits of 
the mass accretion/outflow rates. Based on these conditions 
and using Equation ( |10[ ), we estimate the upper limit of the 
mass accretion rate to be: 


M < 


f Mbh \ . — 1 


8.0 X 10® M0 


M© yr” 


For the canonical values in Equation (111, the mass ac¬ 
cretion rate is M ^ 8 x 10“® M© yr“^ at a distance of 0.4 pc 
from the central engine of NGC 1068. As we only have phys¬ 
ical information of those clouds within the outflow where we 
measure NIR dichroism, the estimated infiow/outflow mass 
accretion rate represents only a fraction of the black hole 
mass accretion rate, M = 0.18 M© yr“^, calculated using 


the bolometric luminosity, \jboi = 9.55 X lO"*"^ erg s ^ (|Woo 


fc Urry|[200^, of NGC10 68. 


Elitzur fc Shlosman (20061 showed that the mass out- 
and mass accretion rate, M, can be 


flow wind rate, M„ 
related as: 
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M, s; 0.14 


T 

1500 K 


Nn 

10^3 cm“^ 




be examined in a sample of AGNs to find general and/or 
extraordinary properties of the outflowing wind model. 


7 CONCLUSIONS 


where Ah is the total column density normalised to 10 ^^ 
cm“^, V is the velocity at the inner radius of the torus nor¬ 
malised to 100 km s“^ based on H 2 O maser VLBI obser¬ 
vations of NGC 1068 at ~0.4 pc by |Greenhill et S!] ( |1996[ ), 
e is the accretion efficiency at the torus radii, and Ii is an 


unknown factor of order unity (Section 6.1 in Nenkova et al. 
||2008b| ). 

We took the to tal c olumn density, Ah = l.dllg'gy x 
10^^ cm~^ (Section 5.1 1 , the velocity of 100 km s“^, the 
accretion efficiency at the torus inner radius of e = 0 . 01 , 
typically used at the torus scales by|Emmering, BlandfordHI 


Shlosman (19921; Elitzur & Shlosman (20061, and we obtain 


an upper-limit of the mass outflow rate of M„ ^ 0.17 Mq 
yr“^ at a distance of 0.4 pc from the central engine of NGC 
1068. 

The torus outflow timescale is tw = Mtorus/Af„. Mtorua 
is the torus mass given by Mtorua = uipAn f N^idV ~ 
10^Ah,23L45U Mq (Elitzur & Shlosman 20061, where L 45 
is the bolometric luminosity normalised to 10 “ erg s“^, and 
Y is the torus radial thickness. We took the bolometric lu¬ 
mino sity to be Lboi = 9.55 X 1 0"^"* erg s“^ (Woo & Urry 
20021, and Y = 6 / 1 ^ (Section 5. l|, and we estimate the mass 


of the torus to be Mtorus = O.TSti'?® x 10^ M©. This value 
is within a factor of three lower than Mtorus = 2.4 x 10® 
Mq using radio observations by Gallimore et al. (1996 1 and 
Mtorua = 2.1(±1.2) X 10® Mq using mm-observations with 
ALMA by Garcfa-Burillo et al. (20141. The difference in 
torus mass can be interpreted as the Clumpy torus htting 
to the IR (1—20 pm) SED is only accounting for those clouds 
emitting in the IR, where a more compact torus would be 
detected. The outflow timescale is estimated to be ^ 10® 
yr. If we assume a torus rotating around the accretion disc 
with a typical Keplerian orbit of Ik = 3 x lO^Mgn 
yr ( Elitzur fc Shlosman||20(i6 1 , the timescale to complete a 
Keplerian orbit at 0.4 pc from the central engine is Ik ~ 10® 
yr. Thus, the obscuring structure of NGC 1068, generated 
by the outflowing wind, can be created in ^100 Keplerian 
orbits, yielding a rotational velocity 5C1228 km s“®. These 
results suggest that the outflowing wind, generated by the 
accretion disc’s magnetic field, can rapidly create an dusty 
obscuring structure around the central engine of NGC 1068. 
Recent time-resolving 3D hydrodynamical models of AGN 
(Schartmann et al 20141 found that the radiation of the 
central engine powers a loss of mass with a rate of 0.1 Mq 
yr“®, with similar inflow mass rates from large scales in or¬ 
der to keep the dusty structure lifetime for a period of ~10 
Myr. In addition, some models assuming gas accretion flows, 


such as adiabatic inflow-outflow solutions (ADIOS, Bland- 


ford & Begelman (19991) and convection-dominated accre¬ 


tion flows (CDAE, Narayan, Igumenshchev & Abramowicz 


( 20001 ), suggest that only a small fraction of the accreted 
matter at the outer radius of the inflow contributes to the 
mass accretion rate at the black hole. This can be explained 
by turbulence and strong mass lost. In this scheme, the ac¬ 
creting efficiency is low (e ~ 0 . 01 ) at the outer radius of 
the inflow, favoring the outflow wind. This scenario should 


We found that the existence and evolution of the torus in 
NGG 1068 can be explained through the hydromagnetic out¬ 
flowing model of AGNs. Assuming Clumpy torus models, 
the polarisation of NGC 1068 in the K' filter is most likely 
arising from radiation of the directly illuminated inner- 
facing clumps of the torus passing through the magnetically 
aligned dust grains located in the low-density outer regions 
of the clumps. We found that the constant component of 
the magnetic field in the plane of the sky is dominant and 
responsible for the dust grain alignment in the torus with a 
strength in the range of 4—139 mG. We presented a direct 
comparison between the P.A. of the NIR polarimetric and 
IR interferometric observations of NGC 1068. Specifically, 
a toroidal geometry is the most likely configuration of the 
magnetic field in the torus. Adopting the estimated magnetic 
field configuration, we And a mass outflow rate of 5:0.17 Mq 
yr“® at a distance of 0.4 pc from the central engine. At this 
rate, the obscuring structure around the central engine can 
be created in a timescale of ^ 10 ® yr with a rotational veloc¬ 
ity of ^1228 km s“®. We conclude that the origin, evolution 
and kinematics of the dusty environment obscuring the cen¬ 
tral engine of NGC 1068 can be explained by a hydromag¬ 
netic outflow wind. Further NIR polarimetric observations 
of a sample of AGNs are required to refine and/or modify 
this approach and find general and/or extraordinary mag¬ 
netic properties in the torus. In addition, MIR polarimetric 
observations are essential to test the effects of dichroism 
at longer wavelengths. Also, mm-polarimetric observations 
with ALMA will allow us to refine intrinsic properties of 
the torus, i.e. dust density, grain sizes, temperature, used 
to estimate the magnetic field strength through the differ¬ 
ent approaches presented in this work, as well as, to test 
the MHD outflow wind model for those clumps at the outer 
edge of the torus. 
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